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^ ■ Abstract 



^ . The problem of bound states in effective field theories is studied. A rescaled 

version of nonrelativistic effective field theory is formulated which makes the 
velocity power counting of operators manifest. Results obtained using the 
rescaled theory are compared with known results from NRQCD. The same 
ideas are then applied to study Yukawa bound states in 1 + 1 and 3 + 1 
dimensions, and to analyze when the Yukawa potential can be replaced by 
a (5-function potential. The implications of these results for the study of 
nucleon-nucleon scattering in chiral perturbation theory is discussed. 
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I. INTRODUCTION 



Effective field theories are an extremely useful tool for studying the dynamics of particles 
at low energies. An effective Lagrangian typically has an expansion in inverse powers of some 
mass scale M, and describes dynamics at momentum scales which are much smaller than 
M. For example, heavy quark effective theory (HQET) [|l]-|^ describes the dynamics of 
hadrons containing a single heavy quark of mass itiq at momentum transfers much smaller 
than TTiQ. The HQET Lagrangian has an expansion in inverse powers of m^, and is used 
to compute hadronic properties as an expansion in Aqcd/''^q, where Aqcd ~ 300 MeV is 
a typical strong interaction scale. The scale at which the HQET Lagrangian ceases to be 
useful is the mass uiq of the heavy quark. HQET can be used to study the interaction of a 
single heavy quark with light degrees of freedom, provided the momentum transfer is small 
compared with mq. 

Systems containing a heavy quark and antiquark (or two heavy quarks) can not be de- 
scribed by HQET. As niq oo, the quark and antiquark form a Coulomb bound state of 
size l/niQasi and typical momentum transfer p ~ asTnq. Perturbation theory is infrared di- 
vergent, with terms of the form (a/w)", where v is the relative velocity of the Q and Q in the 
center of mass frame. These infrared singular terms cause a breakdown of perturbation the- 
ory, and must be resummed. Resummation of the most singular (a/f terms is equivalent to 
solving the Schrodinger equation in a Coulomb potential, and the resummed QQ scattering 
amplitude contains the corresponding bound state poles. One can construct a different ef- 
fective field theory, NRQCD (non-relativistic QCD), that is appropriate for the study of QQ 
bound states in QCD (or its QED analog, NRQED for the study of positronium) The 
terms in the NRQCD Lagrangian are of the same form as those for HQET, but NRQCD has 
a different power counting scheme than HQET. In HQET, the power counting of operators 
is manifest. An operator with coefficient l/mg has a matrix element of order (Aqcd/'^^q)'^- 
Thus, the quark kinetic energy operator (D^/2mQ) Q is of order KQCj^/mQ, and is sub- 
leading. The power counting is more subtle in NRQCD — (D^/2mQ) Q is treated as 
a leading order operator, and higher dimension operators are suppressed not by powers of 
Aqcd/'"^q but by powers of the relative three- velocity of the heavy quarks. We will refer 
to theories such as HQET, in which [iD^) Q is of leading order but (D^/2mQ) Q is 
small, as static theories, and will refer to theories such as NRQCD, in which {iD^) Q and 

(D^/2mQ) Q are both of leading order, as non-relativistic theories. 

In this paper, we study the problem of non-relativistic bound states in effective field 
theories. In Sec. we define a rescaled version of NRQCD (RNRQCD) in which the 
power counting of operators is manifest. Most of the results of RNRQCD follow trivially 
from well-known results in NRQCD, but there are a few advantages, which are discussed 
in Sec. |I[ In Sec. |T| we analyse non-relativistic Yukawa bound states due to the exchange 
of a massive scalar, in both 1 + 1 and 3 + 1 dimensions, and compare the results to those 
in an effective theory in which the scalar is integrated out and replaced by a four-Fermi 
interaction. Naively, the momentum expansion of the four-Fermi effective theory has a 
convergence radius equal to the mass of the scalar, but this is not necessarily the case when 
there are weakly bound states. We study whether the effective theory can "produce its 
own bound state," that is, whether one can obtain a composite weakly bound state in an 
effective theory where the higher dimensional operators have small coefficients. We show 
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that in 1 + 1 dimensions, the effective four-Fermi theory has a radius of convergence of 
order the scalar mass, and produces a weakly bound state with perturbative coefficients for 
higher dimension operators when the parameters of the Yukawa theory are such that there 
is one weakly bound state in the spectrum. If the Yukawa theory has two or more bound 
states of which one is weakly bound, the radius of convergence of the four-Fermi effective 
theory vanishes as the weakly bound state approaches threshold, and the higher dimension 
operators have large coefficients. In 3 + 1 dimensions, the four-Fermi effective theory has a 
finite radius of convergence and higher dimension operators with small coefficients only when 
there is no bound state near threshold in the Yukawa theory. We discuss the implications 
of our analysis for recent attempts to describe non-perturbative aspects of nucleon-nucleon 
scattering by solving the Schrodinger equation for a chiral Lagrangian in which the nucleon- 
nucleon interaction arises from both pions exchange and contact terms [^j-Q. We present 
our conclusions in Sec. TV . 



II. POWER COUNTING IN NON-RELATIVISTIC QCD 

The HQET Lagrangian at leading order in the l/mq expansion is 

£ = Q^ Ud') Q, (2.1) 



where Q is the annihilation field for a non-relativistic quark. It is well known that the 
Lagrangian Eq. (|2.1| ) cannot be used for systems containing two heavy quarks. The basic 
problem arises from the Feynman graph in Fig. |I], when both the intermediate fermions 
are simultaneously almost on-shell. The box graph evaluated in QCD (with propagators 
1/ (^ + m)) has terms which are of order mQ/|k| = 1/ |v|, where k and v are the three- 
momentum and velocity of the external quarks, respectively. The box diagrams in HQET 
can not reproduce this behavior since the Feynman rules are independent of rriQ. As a result, 
QCD can not be matched on to HQET in the QQ sector. The box graph in HQET has a 
loop integral of the form 

d^k 1 1 

(2.2) 



(27r)' k^ + te -k^ + 



The kP integral has a pinch singularity, and is divergent. This singularity is a signal that 
higher dimension operators in HQET are important in the QQ sector. One constructs a 
different effective field theory, NRQCD, with the leading order quark Lagrangian 



£ = gt(,/)0)Q + gt ]g + ... ^ (2.3) 



in which D^/2mQ is considered to be of the same order as . The box graph in NRQCD 
has a loop integral of the form 



d'^k 1 



(27r)^ fco - |k|^ /2mQ + le -fco - |k|^ /2mQ + ie 



(2.4) 



instead of Eq. ( p.2|) . The k integral is now finite, and gives 
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FIG. 1. The box graph in the QQ sector. The double Unes represent nearly on-shell 
non-relativistic quarks and antiquarks, and the curly lines represent gluons. 



— I 



I 



(2.5) 



{2Tlf |k 



which, when the gluon propagators are included, reproduces the mQ/|k| enhancement of 
the box graph in QCD. 

HQET has a simple power counting scheme in which is of order Aqcd- The im- 
portance of various operators is then manifest from the Lagrangian. An operator in the 
Lagrangian of dimension 4 + r has a coefficient of order 1/irtQ, and is of relative order 
(Aqcd/w^q)'^- It is also trivial to count powers of l/mg in loop graphs. The quark and 
gluon propagators are independent of mg, so any Feynman graph with V vertices of order 
l/rrfq . . . ^/itlq has an overall factor of l/rrfq^ where r = J2i=i "^i- NRQCD has a more com- 
plicated power counting scheme, which is discussed in detail in Ref. p. In NRQCD, both 

and D^/2mQ are of the same order, so the quark propagator is i/ {k^ — k^/2mQ + ie). 
The quark propagator depends on mg, so that one can get factors of mg from loop graphs, 
and the power counting rules for loop graphs in NRQCD is not as straightforward as in 
HQET. Furthermore, there are several relevant scales in NRQCD: itlq, the three-moment 
of the heavy quarks rnqv, the kinetic energy of the heavy quarks mqv'^/2, and Aqcd- The 
matrix elements of higher dimension operators in NRQCD are suppressed by powers of v, 
but the V counting is not manifest in the Lagrangian. 

It is advantageous to have an effective field theory with manifest power counting rules. 
One can achieve this for NRQCD by rescaling the fields and coordinates of the usual NRQCD 
Lagrangian. In a nonrelativistic system, E and p are of order rnqv"^ and mqv, respectively. 
Therefore, it is useful to rescale the coordinates so that these are the natural sizes of the 
energy and momentum. Define new coordinates X and T, and new fields and A by 



where A^. = 1/mqv. Requiring 5° and V^/2mg to be both of the same order determines 
At = rnqX"^ = 1/mgf^. This gives the relation between the rescaled energy and momentum 
and K and the original variables and k. 



X = A,X, t = XtT, Q = Ag^, A 



XaoA , A = XaA, 



(2.6) 



= k^rnqv^ 
K = k/mqv. 



(2.7) 



In a nonrelativistic system the rescaled energy and momentum are both of order unity. 




Upon rescaling, the Lagrangian density picks up an overall factor of X^Xt from the change 
of integration variables d^x dt — >■ X^Xt d^X dT, so the (iD^) \1/ term is canonically normal- 
ized if Xq = A~^/^. The gauge field quadratic terms get rescaled to 

^ mgXlXl (V X Af , 
- rriQXlXlo (V^° 

mg \dt J 

In the infrared, A^^ — > oo, and the dominant gauge kinetic terms are (V X A)' and (VA") . 
These terms are properly normalized if 

X^ = X^, = (mQXl)''^\ (2.9) 
With this rescaling the NRQCD Lagrangian becomes 

+^ - vdoAt - gV^fabcA'iA^y 
= vl/t + ^ - -j^A^oT""^ [d^A] - d.A^y + \ {d^A^f + 0{v,gV^) (2.10) 



which will be referred to as the RNRQCD (rescaled NRQCD) Lagrangian. (The effective 
Lagrangian also contains the corresponding terms for the heavy antiquark field.) 

It is clear from the form of the RNRQCD Lagrangian that the effective coupling constant 
in the sector for Coulomb gluons (i.e A^) is a/v, not a. At low velocities, the Coulomb 
gluon interaction must be summed to all orders. Transverse gluons have a coupling constant 
of order g^, and decouple as f — 0. Loop integrals with Coulomb gluons are independent 
of g, V and mg, so the v dependence of a graph may be easily read off from the vertex 
factors. There are no hidden enhancements factors from loop graphs.Q This is in contrast 
with the usual formulation of NRQCD, where 1/v enhancements arise due to factors of 
mg/lkl in the loop graphs. For non-relativistic bound states, v is of order a, and the 
Coulomb interaction in Eq. ( |2.1CI| ) becomes strongly coupled and must be summed to all 
orders. Other interactions, such as those due to transverse gluons, are suppressed by powers 
of V and may be treated using perturbation theory. These well-known results follow simply 
from the scaling of the various terms in the rescaled Lagrangian Eq. (|2.1CI| ). The derivation 
using the original NRQCD Lagrangian is more involved P,p!0 . 



counting for transverse gluons is more complicated, and is discussed in Sec. [I A 
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The effective Lagrangian Eq. ( p.lO|) contains additional higher dimension operators which 



are also suppressed by factors of v. Terms in the effective Lagrangian are relevant, irrelevant 
or marginal, depending on whether the power of v in the coefficient is negative, positive or 
zero.0 Note that in RNRQCD the fields and derivatives are all dimensionless, so terms with 
additional fields or derivatives are not suppressed unless they appear with additional powers 
of V. For example, the operator 

(2.11) 

in NRQCD becomes, in the rescaled theory, 

1 rriQXlXl ^^^4^ ^ ^2 ^y4^^ ^2.12) 

and is of order f^, which agrees with the power counting in Ref. The relation between 
our power counting and that of Ref. is slightly more subtle for operators containing E 
and B fields. For example, the chromomagnetic moment operator 

-^i;^a^-G^,ij (2.13) 

becomes, in the rescaled theory 



-^rnqX^XlXA ^^a^'G^.u'^ = gy^^^a'^'G^.^i/ (2.14) 

whereas in the power counting of Ref. P] this operator is of order v"^. However, the NRQCD 
power counting refers to the size of the matrix element of the operator in a quark- ant iquark 
state. The chromomagnetic gluon must therefore be attached to one of the external quark 
lines, which costs an additional power of g^/v^, hence the matrix element of the operator 
Eq. ( p.l4|) is of order g'^v ~ f ^ in a quarkonium state, as expected. 



A. Tranverse Gluons 

The quark propagator in RNRQCD is 



(2.15) 



/Sfo - |K|V2 + ie 



and the propagator for is 



IK 



2 ' 



(2.16) 



^This is in the renormalization group sense. Irrelevant operators inserted in loop graphs which 
are sufficiently divergent can produce effects that do not vanish as v ^ 0. 
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FIG. 2. Radiative correction to Coulomb scattering. The gluon with a dashed hne is a Coulomb 
gluon. 

which are both independent of v. The transverse gluon propagator is 

- ^ { - ^ 2 • (2-17) 

V ' |K|' I |K|' - (i^o)2 V ) 

which depends on w, so there are potential 1/v enhancements from graphs with internal 
transverse gluons. A generic loop integral can be evaluated by first doing the integral 
using residues. The residue of the transverse gluon propagator at the pole |K|/t> is 
— 1/f |K|, which is enhanced by 1/v. The transverse gluon propagator contribution from 
other poles (such as fermion poles) is typically of order unity, since these poles are at values 
of of order unity. Thus transverse gluon loops can have 1/v enhancements from regions 
in the momentum integral where the transverse gluon is on-shell. Consider, for example, 
the graph in Fig. |^, where the gluon in the loop is a transverse gluon. The diagram has a 
1/v enhancement from the region in loop integral where the fermion and virtual gluon are 
on-shell, i.e. from physical tranverse gluon radiation. 

The transverse gluon coupling constant is g\/v, which is a factor of v smaller than that for 
Coulomb gluons. Thus transverse gluon loops obey naive v counting, and are v"^ suppressed, 
unless the cut part of the graph contributes to transverse gluon radiation. In the latter case, 
the transverse gluon graph has a 1/f enhancement over the naive w-counting, and is only 
suppressed by one power of v. 

It is not possible to choose a rescaling scheme which has a manifest f-counting scheme 
for real and virtual transverse gluons. The reason is that quarks behave like non-relativistic 
particles when transverse gluons are exchanged between them, but like static particles when 
one of them radiates an on-shell transverse gluon. 



B. Coulomb Scattering 

At leading order in w , the only diagrams which contribute to QQ scattering in the effective 
theory are the Coulomb ladder graphs of Fig. ^ where the propagators are denoted by 
gluons with dashed lines. Crossed ladder graphs such as Fig. ^ vanish when both gluons are 
Coulomb gluons. The one-loop box graph has an integral of the form 

r dK° d^K 1 1 1 x (2 18) 

J (27r)(27r)' IK^K + P/ - Pi|' (Pi - 7^)° - |Pi - K| V2 + «e 
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FIG. 3. The leading order contribution to QQ scattering. The gluons with dashed hnes repre- 
sent Coulomb gluons. 




FIG. 4. The crossed box graph. 



(P2 + i^) - |P2 + K|72 + ie 
The integral can be done by contour integration to give 

dK^ 1 1 



(27r) (Pi - Kf - |Pi - K|72 + {P^ + Kf - + K]' /2 + le 

i 



(2.19) 



(Pi + P2)° - |Pi - K|' /2 - IP2 + K|' /2 



which is the Schrodinger Green function of non-relativistic quantum mechanics. The ladder 
graphs of Fig. ^ can be evaluated by first doing the K'^ loop integrals. It is easy to see that 
each loop gives the Schrodinger Green function for the intermediate two-fermion state. The 
sum of all the ladder graphs gives the non-relativistic Schrodinger equation in momentum 
space for a Coulomb potential. This result is, of course, well-known. What is different here is 
that in the rescaled theory, the sum of the leading order graphs is identical to the Schrodinger 
equation. While this is also the case for NRQCD in Coulomb gauge, the rescaling allows 
similar results to be obtained in non-gauge theories, such as Yukawa theory (which we will 
discuss at length in the next section), where there is no gauge freedom in choosing the 
propagator. 



C. Heavy-Light Systems 

One can see explicitly why D^/2mQ is a relevant operator for heavy-heavy bound states, 
but not for heavy-light bound states. Consider a non-relativistic bound state of two particles 
with different masses tuh and m^,, with tuh 3> m^. Applying the rescaling Eq. ( |2.9|) with 
mq rriL, one finds that the fermion kinetic terms are 

V^Iz^Vl + ^V^IV^^l (2.20) 
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for the light fermion, and 



2m 



H 



(2.21) 



for the heavy fermion. The ip'^'W'^ip operator is comparable to ip'^d^ip for the light particle, 
but is smaller by niL/mH for the heavier particle. The heavier particle can be treated as a 
static source (as in HQET), but the lighter particle must be treated as non-relativistic (as 
in NRQCD). 



D. Scaling Dimensions 

The velocity scaling rules have a renormalization group interpretation. In the non- 
relativistic limit, the scaling dimensions of space and time should be chosen to be 



[x] = -1, [t] = -2, 
instead of the usual choice in relativistic theories 

[x] = -1, [t] = -1. 
With the choice Eq. (g^ , one finds that 



[C] = 5, 



3/2, 



A' 



3/2, [A] = 3/2, 



(2.22) 



(2.23) 



(2.24) 
(2.25) 



fcO] = 2, \k] = 1. 



Operators are relevant, irrelevant, or marginal depending on whether their dimensions are 
less than, greater than, or equal to [C] = 5. For example, the Coulomb interaction term 
\|/t^Ovl> has dimension 9/2 = 5 — 1/2, so the Coulomb interaction is relevant, and of order 
y-i/2_ rjj_^g transverse gluon interaction term is \l/^p ■ A\I' and has dimension 11/2 = 5 + 1/2, 
so the term is irrelevant, and of order f This agrees with the powers of v in the Lagrangian 
Eq. (CTD - 



III. POWER COUNTING IN CHIRAL PERTURBATION THEORY 

There has been much recent interest in applying the techniques of non-relativistic effective 
theories to nucleon-nucleon scattering 0-0. The goal is to describe low-energy nucleon- 
nucleon scattering using a chiral Lagrangian where massive excitations (such as p's and cj's) 
have been integrated out and replaced by an effective four point interaction. Nonperturbative 
effects, such as the large scattering length in the ^5*0 channel or the deuteron in the ^Si 
channel, could be described by solving the Schrodinger equation in the effective theory. 

The rescaling of the previous section can be extended to chiral perturbation theory for 



heavy nucleons The leading terms in the nucleon-pion chiral Lagrangian have the 

form 
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^-fW-^W (3.1) 

2 ^1 /-'t./A'-^r./.t.- a2 



describing the interaction of a pseudoscalar Goldstone boson cf) with decay constant and 
couphng (jf, where we have suppressed flavour indices. The leading effects of massive excita- 
tions are contained in the dimension six operators {ip'^tl))'^ and [ip'^aipY with dimensionless 
coefficients Ci and C2. In 3 + 1 dimensions the appropriate rescaling is 

X = A^X, t = MXlT, ip = A;^/2^, = (mA^)~^^^ $ (3.2) 
giving the rescaled Lagrangian 

£ = vl/t (^^^' + ^) * - V-l- ■ * Vv& (3.3) 

/<^n^\2 1/^^x2 CiM^i; /^x \2 021/^1' /^x 

- 2 ^^"^^ ^ (^^) ^ ) • 

Both the derivative interaction and the dimension six operators are irrelevant in the infrared. 
Each pion exchange in a Feynman diagram contributes a factor of g'^M^v/ f^, which has the 
same dependence of M and v as an insertion of one of the dimension six four-nucleon 
operators. This reproduces the power counting of Ref. 0: a graph with n insertions of 
dimension six operators and m ladder pion exchanges scales as 



(3.4) 



where Q = Mv is the three momentum of the scattering nucleons (in the centre of mass 
frame). Furthermore, with this rescaling the power counting of Ref. 0, in which poles in 

in the pion propagator were treated as higher order in the nonrelativistic expansion, is 
manifest. Since pion emission is kinematically forbidden as f ^ 0, the temporal piece of the 
pion kinetic term may be treated as a higher-order insertion. 

Since the pion-nucleon interaction terms are irrelevant in the infrared, unlike the interac- 
tions in QED and QCD, perturbation theory for nucleon-nucleon scattering does not break 
down at threshold for weak coupling. Multi-loop bubble graphs are not enhanced by powers 
of 1/v, but instead suppressed by powers of v. It was argued in Refs. that the large 
scattering length in the ^5*0 channel and the deuteron in the ^Si channel in nucleon-nucleon 
scattering are signs of the breakdown of perturbation theory, due to the fact that for large 
M the factor in Eq. ( ^^ is not small. It was proposed by these authors that the appropriate 
description of nucleon-nucleon scattering could be obtained by summing all terms of order 
(QM)"; corrections to this would be suppressed by powers of Q/M, and could be calculated 
systematically. 

However, as noted in Ref. , the validity of this approach rests on the assumption that 
coefficients of higher dimension four-nucleon operators in the chiral Lagrangian with n spatial 
derivatives are smaller than M". Otherwise, the effective theory description would break 
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down, since the contributions from higher dimension operators with unknown coefficients 
would be as large as the graphs which are being summed in the effective theory. In it was 
argued that this assumption was consistent in the ^5*0 (but not in the ^Si) channel, since 
no divergences were found in the Feynman diagrams contributing to ^5*0 channel scattering 
which would require counterterms scaling like (QM)". 

In this section we will investigate the validity of this power counting in a simplified model. 
We will neglect the Goldstone bosons, and simply consider a theory of a non-relativistic 
fermion of mass M coupled to a scalar of mass m <C M, 



Y 



('^° ^jPl^ - ^'/'V^V + I (5V)' - I (V0)' - Im'cp'. (3.5) 



The scalar plays the role of the p, uj and other excitations. At low momenta, p <^ m, the 
scalar may be integrated out, resulting in an effective four-fermi theory 



2M 



j + + (3.6) 



where the ellipses represent higher dimension operators, suppressed by powers of p/m. The 
Yukawa potential arising from scalar exchange in Eq. (|3.5|) is replaced by a series of delta 
functions and their derivatives in the effective theory. We will refer to these two theories as 
NRY and NR(5 for non-relativistic Yukawa, and non-relativistic 5-function, respectively. 

We now wish to ask whether NR.6 correctly describes the V'V' scattering amplitude of 
NRY at low momentum p < Pmax, in a parameter regime where the Born approximation 
fails. By low momentum, we mean p < Pmax fn, where Pmax is held fixed as one varies the 
parameters in the Yukawa theory. This is analogous to the question of whether the dimension 
six operators in Eq. (|3.3|) correctly describe nucleon-nucleon scattering at momenta much 
smaller than rrip, in a regime where the Born approximation fails. We consider the problem 
in both 1 + 1 and 3 + 1 dimensions. In 1+1 dimensions, we will show that NKS correctly 
describes the scattering for weak coupling; however, for a coupling large enough that more 
than one bound state exists the higher dimension operators in NKS become important, and 
the effective theory breaks down below Pmax- In this case, the full theory is required to 
correctly describe the low-energy physics. In 3+1 dimensions, there are no bound states 
for sufficiently weak coupling. In this case NR5 correctly describes the scattering and no 
resummation of ladder graphs is necessary. However, when the coupling is strong enough 
to form bound states, the situation is analogous to the 1+1 dimensional case with excited 
states. Once again the higher dimension operators become important, and again the effective 
theory description breaks down below Pmax- We find by explicit calculation that there are, in 
fact, higher dimension operators in the low energy theory which scale like (QM)", and spoil 
the power counting of Refs. f^-^. We will comment on the relation between our results and 
the effective range expansion, as well as the case where Pmax is allowed to vary, in Sec. |111 C| . 

A. 1 + 1 Dimension 

In 1 + 1 dimensions, the mass dimensions of the fields and couplings in Eq. ( |3.5| ) are 

[^] = 1/2, [0] = 0, [M] = 1, [m] = 1, [g] = 1. (3.7) 
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The Yukawa potential due to scalar exchange is 



= / tI^^^ = -|^e-N. (3.8) 

An attractive potential in 1 + 1 dimensions always has a bound state. The Schrodinger 
equation for a Yukawa potential is not analytically soluble. However, when g is small, the 
state is weakly bound and spread out over a large region. In this limit the Yukawa potential 
can be approximated by the 5-function, 

_ -m\x\ _^ -^S(x) (3.9) 
The 5-function potential has a bound state with energy E = —g^M/Am^, and wavefunction 



g^M 



ij{x) = v^e-"l"l, K = (3.10) 

Corrections to the (5-function result will be suppressed by powers of K/m, the ratio of the 
size of the state to the range of the potential. For /t/m ~ 1, excited bound states appear, 
and the 5-function approximation completely breaks down. 

In the effective field theory language, this indicates that NR5 will only correctly describe 
bound states for weak coupling, n/m <^ 1, and will break down at larger values of the 
coupling. Since an excited state may be arbitrarily close to threshold, the effective theory 
may therefore break down for tpip scattering at arbitrarily low energies, which is not the usual 
behavior one expects from a low-energy effective theory, where the importance of higher 
dimension operators is set by powers of p/m. In this case, the effects of higher dimension 
operators must be suppressed instead by powers of n/m. Therefore, it is instructive to 
analyze this problem by matching NRY onto NR5 and seeing how the scale k enters the 
problem. 

To make the power counting manifest we use a rescaled Lagrangian as in Sect. |l[ In 
1+1 dimensions, the appropriate rescaling is 

X = A^X, t = MXlT, ^ = A;^/^^, (j) = (MA^)"^/^*. (3.11) 

where A^ = 1/Mv, giving 

^ I 2 M 2 ^ ^ 2^ ' 2A'P ^ ^ 



The powers of v can also be obtained using modified scaling rules, as in Sec. |II D| . The 
modified scaling dimensions are 



[x] = -1, [t] = -2, [C] = 3, [ij] = 1/2, [</.] = 1/2. (3.13) 



[C] = —3/2 and [0^] — [C] = —2, so that these terms scale as v '^/^ and 



V' 



-2 



respectively. Since the temporal piece of the scalar kinetic term in Eq. ( |3.12| ) is of order 
f^, it may be neglected at the order at which we are working. Thus RNRY (rescaled NRY) 
describes a fermion coupled to a static scalar field, with coupling constant gv~^^'^/M. The $ 
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interactions must be summed to all orders, since the Yukawa coupling is large. Because the $ 
propagator has no poles in (emission of scalars is kinematically forbidden as f — 0) only 
ladder graphs are non- vanishing in the nonrelativistic theory. For simplicity we consider the 
scattering of two off-shell fermions with energy E and zero momentum, since the scattering 
amplitude will have the same singularities as for on-shell states. We will also treat the V^'s 
as distinguishable, since the additional graphs from Fermi statistics are not relevant to our 
arguments. The tree-level $ exchange graph in the ladder sum in Fig. ^ gives 

,2 
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1 



,2 



9 



(3.14) 



where jjr = m /M v . The box graph gives 



Y 



9' 



4 



dK 



9 



2/i + 



4/i3yi(^ + ^)^ 
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(3.15) 



where —e = {Pi + P2)° is the total rescaled energy of the incoming particles. The two-loop 
box graph gives 



rY 



M^v^ J {2n 



dK r dL 



1 



1 



(27r) (7^2+^2) (L2 + ^2)((^_^y 

1 



(Pi + p^y - K^ (Pi + P2r - 

6/1^ + 26 



'M^v^ 12/i% (/i + (^ + 2^) (2/i + ^) 



■ 9 



— I 



(3.16) 



The 0{l/v) term from the box graph Eq. (|3.15|) is down by a factor oi n/m relative to 
the 0{l/v) term due to tree-level exchange, and so may be neglected in the weak coupling 
limit, K <^ m. In this limit, one may sum the most singular contributions of the graphs in 
Fig. ^ to obtain for the ladder sum 



lA 



rn?v 



1 + h 

2ni?v 



2'im?v 



+ ... 



7ii?v 1 — g'^e ^1"^ I {2vn?'v) 



(3.17) 



which has a pole at e = g^ j (4m^f^). Rescaling back to physical units gives E = 
— (y'^M/(4m'^), which is the correct bound state energy for n/m ^ 1. There is always 
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+ 



FIG. 5. Ladder graphs in the Yukawa theory. 



at least one bound state pole, even for weak coupling, because the box graph diverges as 
V ^0. 

When K/m is not small, the ladder sum is no longer a geometric series, since an n-loop 
ladder graph has singularities of the form 



K Y 1 

mj v' 



n-l 



1 1 



V V\ e 



(3.18) 



(ignoring factors of g/m) as can easily be verified by a dimensional estimate of the n-loop 
graph. This behavior can also be seen from the explicit computations Eqs. (|3. 14 - 3. 16 ). Each 
term in Eq. ( p.l8|) is as large as a term in the sum Eq. (|3.17|) , and can not be neglected. 
The scattering amplitude will be given by the Green function for a Yukawa potential (which 
is not known analytically) and will have poles at the energies of all the bound states. 

We now wish to integrate the $ field out of the theory to obtain the rescaled effective 
four-Fermi theory]^ 



+ .... 



The tree level matching condition arises from $ exchange and gives 



(3.19) 



(3.20) 



where we denote by /i*^^'-* the r-loop contribution to h. The one-loop matching condition is 
shown in Fig. ^ The box graph has already been evaluated in the full theory in Eq. ( p. 151) . 
In the effective theory, the box graph is 



dK'^ dK 



1 



1 



(27r) (27r) (P^ - Kf - + %e {P^ + Kf - + 



le 



dK 



(27r) K^ + e 

"1/2 



(3.21) 



Using the matching condition Eq. (|3.2CI|) , we see that this reproduces the leading term in 
Eq. ( [3.151 ). It must do this, because one can not write down an interaction proportional to 



The V dependence of 

(^$t$^ 

can also be obtained using Eq. (|3.13D , since 



-1. 
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an inverse fractional power of e in the low energy RNR(5 effective Lagrangian. The 0{l/v) 
term in Eq. ( |3.15| ) is reproduced in the effective theory by the one-loop matching condition 



to /i, 

h^^^ = --^^ = h^^\ 3.22 

4 4 m 

At two loops, one needs to compute the matching conditions in Fig. ^. The two-loop graph 
in the effective theory with two insertions of h^^^ reproduces the 1/e term in Eq. ( p. 161) , 
and the one-loop graph with one insertion of h^'^^ and one insertion of h^^^ reproduces the 
1 / ^/e term. The matching correction to h at two loops requires knowing the coefficients of 

operators of the form (^\E'"''V\1'^ . This is a complication of the 1-1-1 dimensional analysis 
which is not present in 3 -|- 1 dimensions. 

The tree- level matching condition for h only dominates when /t/m <^ 1. In this limit, 
the sum of bubble graphs in NR(5 shown in Fig. |^ again gives a geometric series, 

1 

which reproduces the result Eq. ( |3.15| ) of RNRY for k <^ m. 

When K/m is not small, we have already seen in the full theory that the scattering 
amplitude is not given by the simple geometric series ( |3.17| ); the same is of course true in 
the effective theory. Consider a four-fermi operator in the effective theory containing one 
time derivative, 




FIG. 8. Two loop matching condition from NRY to NR(5. 
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—iP^OqiIj^Ij^^ + h.c. (3.24) 
m 

Like h, higher loop matching conditions to c will be suppressed by powers of /t/m, so the 
matching will have to be done to all orders when k/iti is not small. In the rescaled theory, 
this operator becomes 

l^^t^P^^t^ + h.c. (3.25) 
V m 

In the Born approximation the matrix element of this operator may be neglected compared 
to the terms in the series ( 3.17| ); however, when it is inserted into a two-loop bubble diagram. 



it gives a contribution to the scattering amplitude proportional to 

c Mf ^ / K \ 2 c 



f 2 t> ra \m 



(3.26) 



which, for n/m ~ 1, is the same size as the first term in the geometric series. Similarly, 
operators of arbitrarily high order will contribute at all orders in 1/v to the scattering, 
suppressed only by powers of k/ui, so the effective field theory is unable to describe bound 
states in the theory when this parameter is not small. An infinite number of operators is 
required in the effective theory to reproduce the scattering amplitude in the full theory. 
A simple way to see that the effective theory breaks down is to note that in RNRY, the 
effective mass of the $ is m/\/2Mv ~ m^/g'^M = m/n, for v ~ g^/m?. Thus, for k <ti m 
the $ field is heavy, and may be integrated out, with higher dimension operators suppressed 
by powers of Mv/m ~ n/m. For k m, the $ is light, and cannot be integrated out of the 
theory if a light bound states is to be properly described. 

B. 3 + 1 Dimensions 

The analysis of the previous subsection can be repeated in 3 + 1 dimensions. We will use 
the same symbols as in 1 + 1 dimensions, but they now have dimension 

[tfj] = 3/2, [0] = 1 [M] = 1, [m] = 1, [g] = 0. (3.27) 

The three dimensional Yukawa potential is 

^ J (27r)^ |q|' + m2 47r |x| ' ^ ' ' 

The Schrodinger equation with a three dimensional Yukawa potential does not necessarily 
have a bound state. The bound state first appears when 

^^>1.7 . (3.29) 
47r m 

Under the rescaling (|3.2| ), the NRY Lagrangian in 3 + 1 dimensions becomes 

£ = vl/t (",50 ^ ^ _ 4=$vl>tvl> + ^ (5°$)' - - {V^f - -^^$2 (3.30) 
V 2 J y/v 2 V y 2 ^ ^ 2 M^v^ ^ ^ 
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while the rescaled version of NR5 gives the RNR(5 Lagrangian 



The modified scahng rules of Sec. [II D| give 

[x] = -1, [t] = -2, [C] = 5, 



3/2, 



3/2, 



(3.31) 



(3.32) 



so that 



^ip^ip] - [C] = -1/2, and (t^V 



- [£] = 1, which gives the v dependence of the 

coefficients in Eqs. (|3.30 , 3.31 ). In three dimensions, the Yukawa interaction is a relevant 
operator, and the four-Fermi interaction is an irrelevant operator. This is a sign that the 
four- Fermi theory will have trouble correctly reproducing the behavior of the Yukawa theory 
when there is a light bound state in the spectrum. 

Tree-level $ exchange in RNRY gives the amplitude for \['\E' scattering 



and so the tree level matching condition is the same as in 1 + 1 dimensions, 

^2 



(3.33) 



(3.34) 



The one-loop box graph in RNRY is 

^4 " d^K 



Y 



■ 9 



1 



(27r)-' |Kr + |K + Pi - P'll^ + /i2 
1 



X 



« 9 



{Pi + P2) 
4 1 



Kr/2- |Pi + Kr/2 



(3.35) 



m 



for off-shell incident particles with Pi = P2 = 0, where as before e 
one-loop bubble graph in the low-energy theory is 



(Pi + P2)°. The 



1 



(27r) (27r)^ (Pi - K 
1 



'\0 



|Pi -K|V2 + ie 



P2 + Ky -\P2 + K\' /2 + ie 
d^K 1 



{2Tiy\Ky + e 

,1/2 



(3.36) 
(3.37) 
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where we have evaluated the integral in dimensional regularization. Note that an integral 
which diverges as K"^ , with r an odd integer, is finite when regulated with dimensional 
regularization, and needs no subtractions. The result Eq. ( p.37| ) correctly reproduces the 
leading non-analytic term in Eq. ( p.35| ). However, the analytic term is absent, so one needs 
to add a one-loop matching contribution to 



g^M _ g'^M 1 
Svrm^ 47rm 2 



(3.38) 



In the region of parameter space for which bound states exist in the Yukawa theory, this 
term is at least as large as the tree-level matching term h'^^^ , and so all orders in the loop ex- 
pansion must be calculated in order to calculate the matching for h. The complete matching 
condition for h will be proportional to the scattering length for a Yukawa potential, which 
is not analytically known in closed form. Furthermore, there is no reason for h to be of 
the naive size 1/m?. In particular, if there is a bound state near threshold, h will be much 
larger. 

The two-loop ladder graph in the full theory is 



(47r)2T;3 (^2 



In 



e] 



(2/i + eV^)^ 
3/i(/i + 2£V2) 



. g^M^v 



3 



(47r)2m4 

the two-loop bubble graph in RNR5 is 



Am? 



7 + 8 In ■ 



(3.39) 



3 ^6 ^3 



(47r)2m^ 

and the one-loop graph with a single insertion of h^^^ is 



(3.40) 



4„,2 



a/2 



9' 



1/2 



An (47r)2m4 
The sum of the graphs on the right hand side of Fig. ^ is 

. g^M^v 



m 



(3.41) 



+ . „ 4 



m 



Aw? 



(3.42) 



Comparing Eq. ( p.42|) with Eq. (|3.39| ), we see that the nonanalytic term of order e^i"^ is 
again reproduced in the effective theory, as it must be, but that the terms of order e and £° 
terms are not reproduced. Therefore there is a two loop contribution to 



(47r)2m4 3 \Anm^ 



3' 



(3.43) 



which once again is at least as large as h^'^^ in the region of interest, while at order e the 
difference between the graphs in the two theories 
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,^£f^f3 + 8lnil (3.44) 



{AnYm^ Am? \ 3 
contributes to the matching conditions of an operator such as ^Hd^^^^'^ , with coefficient 



f— — — - + 21n- = - + 21n- h^^^M^v. 3.45 

(47r)2m6 V4 3/ V4 3/ V ^vr / ^ ' 

The important observation is that this counterterm contributes to the e term at order 
(QM)^, the same order as the two-loop bubble graph in the effective theory; hence, without 
knowing the coefficients of the dimension eight operators in the effective theory all terms of 
the form {QM)^ are not summed. It is clear that a similar situation exists at higher loops: 
each loop graph in the effective theory is the same order as a counterterm, and so without 
knowing the counterterms, the graphs may not be summed. Thus, without including all 
higher dimension operators, the effective theory does not correctly sum all terms of order 
(QM)". 



C. Bound States and The Effective Range Expansion 

Let us look for bound state poles in the scattering amplitude. In 1 + 1 dimension, the 
bubble chain sum in the effective theory had contributions proportional to inverse powers 
of V, which diverged in the low energy limit. Thus, for weak coupling there was a bound 
state at threshold, which was reproduced in the low-energy theory. In 3+1 dimensions, the 
full theory with a Yukawa potential has a bound state at threshold when g'^M/Anm ~ 1.7. 
Since in the bubble sum in NR(5 only positive powers of v occur, all higher order graphs 
contributing to the scattering amplitude vanish at low energy. Therefore, there can be no 
bound state pole at threshold unless one of the coefficients in the effective Lagrangian NR(5 
diverges. In particular, there can be a bound state at threshold if /i ^ oo. To see this in 
the 5-function effective theory, one must evaluate the matching condition for h to all orders, 
since g'^M/Anm is not small. 

For finite h the bubble sum in the effective theory may be summed to find a bound state 
for finite h. Neglecting the contributions from operators of dimension > 6, the bubble chain 
sum is 



ihM^v 



Ait \ A-n: 



ihM'^v , , J . (3.46) 

1 + hMV/^/AiT ^ ' 



which is the expression given by Weinberg . There is a pole at 



All 

or rescaling back to physical units, at 



1, (3.47) 
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Air 



-1. 



(3.48) 



which, for < 0, has a bound state below threshold. The bound state energy is 

E-. 



(3.49) 



When \h\ oo, the bound state approaches threshold. 

We have argued, however, that the bubble sum in the effective theory does not correctly 
sum the terms of order [QM)^. What relation does this bound state then have to a bound 
state in the full theory? To answer this, it is useful to recall some results of the effective 
range expansion for potential scattering. 

The scattering amplitude iA is related to the phase-shift 6 by 



k cot 6 



Ait 



+ ik 



(3.50) 



In the effective range expansion, k cot 6 is expanded in powers of k, 

kcot6 = -^ + ^rek^ + C4k^ + .... (3.51) 
This expansion gives us some useful additional information, because it is known to have a 



radius of convergence > m/2 [|T2|,[T3|. Here a is the scattering length, and re is the effective 
range. The first two terms in this series provide a very good approximation to the measured 
nucleon scattering cross-sections. As shown in Ref. 0, Eqs. ( |3.50| ) and ( |3.51| ) imply that 
the four-Fermi operators in NR(5 have a momentum expansion of the form 



Air 
M 



a H — a rf,k + a -r^ — a Ca] k + 
2 V 4 



(3.52) 



Rescaling these coefficients, we find that h is related to the exact scattering length a of the 
Yukawa theory. 



An 

h = — —a 
M 

whereas a two-derivative term like \]/tV^\&\E'^\E' has a coefficient of order 



(3.53) 



(3.54) 



The scattering length diverges when there is a bound state at threshold, so that h — »• 
oo. It is known from potential scattering that the other coefficients in the effective range 
expansion (|3.51| ) do not diverge |1I2|JT3[| . In this case, the coefficients of the higher dimension 
operators diverge as well, as is evident in Eq. (|3.52|), and the divergent behavior of the various 



coefficients is highly correlated. It is difficult to see this behaviour by studying the matching 
conditions to the effective theory, since one must perform the matching to all orders in the 
loop expansion. However, the effective range expansion gives nontrivial information on the 
form of the complete matching conditions. 
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If we now consider ipip scattering, the two-bubble graph from the dimension six operator 
gives a term of order 



4:71 



3 

(3.55) 



which is to be compared with the counterterm contribution ( |3.54| ). The terms have the same 
exphcit dependence on M and v, as we have aheady argued, so the bubble graph is the same 
order as the counterterm using the power counting scheme of P-^. Furthermore both terms 
are suppressed by two powers of Q = Mv relative to the contact term, so even when M is 
large it is not necessary to resum the bubble graphs in this scheme. However, the bubble 
graph will dominate when 

> 1. (3.56) 



Summing bubble graphs in the effective range expansion therefore corresponds to an expan- 
sion in powers of a/47rre, not Q/M. Therefore, when the scattering length is large (and 
therefore the bound state is nearly at threshold) the sum of the bubble graphs in NR5 
correctly describes the bound state, and ( |3.49| ) is valid, up to corrections of order inre/a. 
However, as discussed in 0, when the scattering length is large the effective theory breaks 
down not at A; ~ m ~ r~^, but at a much lower scale 




(3.57) 

which, in the region where the bubble graphs dominate, is much less than m. The region of 
convergence is not controlled by the scalar mass, and vanishes as the bound state approaches 
threshold. Thus, we still do not have the situation we had in weakly bound 1+1 dimensional 
theory in which the bound state was correctly predicted in the low-energy theory but the 
radius of convergence of the theory was set by the heavy scalar mass. 



D. Implications for Chiral Perturbation Theory 



The results are unchanged when pions are included in the chiral Lagrangian, since pion 
exchange scales with M and v in the same way as the dimension six contact interactions. 
We may therefore conclude from this simple model that chiral perturbation theory cannot 
sum all terms of order (QM)" without including an infinite number of higher dimension 
operators. Using the effective range expansion, we see that the bubble graphs only correctly 
describe the low-energy bound state when the effective field theory breaks down at a scale 
much less than the symmetry-breaking scale A^sb (the analogue of the scalar mass in our 
simple model). Therefore the standard form of heavy nucleon chiral perturbation theory 
does not provide the correct description of nucleon-nucleon scattering up to energies set 
by the scale of the heavy excitations which have been integrated out of the theory. This 
does not mean that the effective field theory idea is useless: it does mean, however, that the 
usual chiral Lagrangian is not the correct effective field theory. If the appropriate low energy 
degrees of freedom are introduced by hand as new fields in the effective Lagrangian, it should 
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be possible to correctly describe low energy nucleon-nucleon scattering in an effective field 
tlieory.0 However, the properties of the low energy degrees of freedom are not determined 
by the parameters in the standard chiral Lagrangian, unless all higher dimension operators 
are included. 

IV. CONCLUSIONS 

Rescaled non-relativistic effective field theories simplify the study of weakly bound states 
in quantum field theory. The use of RNRQCD makes the v power counting of NRQCD 
manifest, and reorganizes the perturbation expansion in a more systematic way. Rescaled 
effective theories can also be used to study Yukawa bound states due to scalar exchange. 
In 1 + 1 dimensions, a bound state occurs at weak coupling because the n-loop ladder 
graph diverge is proportional to l/(— i?)"/^. In this limit one may sum the diagrams of 
Fig. and obtain the properties of the bound state, without worrying about higher loop 
matching conditions. In three dimensions, one needs a critical coupling before there exists 
a bound state, and the problem is intrinsically strongly coupled. One needs to include 
the full Yukawa interaction to study the bound state. Equivalently, the matching from 
NRY to NR5 must be performed to all orders in the loop expansion, and the higher order 
matching terms are relevant for the bound state. We have shown by explicit computation 
that the power counting scheme of Refs. fp-^] that involves summing powers of QM does 
not hold for Yukawa theory. One can not, in general replace the Yukawa potential by a 
5-function potential to study even weakly coupled bound states in three dimensions, except 
for a very limited region in momentum space |k| < y^2/are, which vanishes as the bound 
state approaches threshold. 
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